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Abstract 

We study the Littlest Higgs model with T-parity at the LHC through pair 
productions of the T-odd top quark partner (T_) which decays into the top 
quark and the lightest T-odd particle. We identify the region of parameters 
favored by the electroweak and cosmological considerations. The signal and 
background events are simulated with fast detector simulation to study the 
discovery potential at the LHC. We find that the hemisphere analysis recently 
proposed by the CMS collaboration is very useful to separate the signal from 
the ti background. We discuss the observability of the top tagged signal in 
the effective mass (M e £j) versus the transverse missing energy (-Ermiss) plane. 
We show that, for all our sample parameter sets with Mt_ < 900 GeV, the 
excess of the signal over the background can be visible as a bump structure in 
the -E/rmiss distribution for 50 fb _1 at relatively high M e g intervals. 
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1 Introduction 



The Higgs sector in the Standard Model (SM) receives large quadratic mass correc- 
tions from top and gauge boson loop diagrams. New symmetries involving top-Higgs 
and gauge- Higgs sectors below 0(1 TeV) are proposed to remove the corrections. 
One of the physics targets of the ATLAS [1] and CMS [2] experiments at the LHC 
is to find new particles predicted by such symmetries. 

The most important new physics candidate is the Minimal Supersymmetric Stan- 
dard Model (MSSM)|3j. This model predicts quark and lepton partners with spin 
(squark and slepton) and those of gauge and Higgs bosons with spin 1/2 (gauginos 
and Higgsinos). Thanks to the cancellation among boson and fermion loop diagrams, 
the quadratic corrections to the Higgs mass term completely vanish. At the LHC, 
strongly interacting supersymmetric particles such as squark and gluino will be co- 
piously produced. They decay into relatively light electroweak (EW) superparters, 
such as the chargino, neutralino, and slepton. The lightest supersymmetric particle 
(LSP) is stable due to the R-parity of the model. The decay products of the squark 
and gluino contain at least one LSP. The LSP escapes from detector without any 
energy deposit, giving a missing momentum signature to the events. The signature 
of supersymmetric particles has been studied intensively by many groups [HOH]. 

Alternative scenarios which do not rely on supersymmetry to cancel the quadratic 
divergent corrections have been proposed and studied. The little Higgs model [5] is 
one of these alternatives. In this model, the Higgs boson is regarded as a pseudo 
Nambu-Goldstone boson, which originates from the spontaneous breaking of a global 
symmetry at certain high scale, and the global symmetry protects the Higgs mass 
from the quadratic radiative corrections. The simplest version of the model is called 
the littlest Higgs model |6]. The global symmetry of the model is SU(5), which is 
spontaneously broken into SO (5). The part of the SU(5) symmetry is gauged, and 
the gauge symmetry is SU(2)i x SU(2) 2 x U(l)i x U(l)2. The top sector is also 
extended to respect the part of the global symmetry. 

This model, however, predicts a large correction to the E W observables because 
of the direct mixing between heavy and light gauge bosons after the EW symmetry 
breaking. The precision EW measurements force the masses of heavy gauge bosons 
and top partners to be O (10 TeV), reintroducing the fine tuning problem to the 
Higgs mass [7j. A solution of the problem is the introduction of T-parity to the 
model which forbids the mixing j8]. This is the symmetry under the transformations, 
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SU(2)i «-> SU(2) 2 and U(l)i «-> U(l)2- All heavy gauge bosons are assigned a T-odd 
charge, while the SM particles have a T-even charge. The matter sector should be 
extended so that T-odd partners are predicted. The lightest T-odd particle is the 
heavy photon, which is stable and becomes a candidate for dark matter [5]- [11]. 

Starting from a different underlying theory, the littlest Higgs model with T-parity 
ends up predicting a similar phenomenology to that of the MSSM. We view that this 
is indispensable for the model to remove the hierarchy problem. One needs a new 
symmetry to protect the Higgs mass to reduce the quadratic divergence of the theory 
in a meaningful manner. The symmetry must involve both the top quark and the 
gauge sectors, because the Higgs couplings to these particles are the dominant source 
of the divergences. Unless some parity is not assigned to the gauge partner, large 
corrections to the EW observables are expected, and those are not acceptable after 
the LEP era. Some of the top and gauge partners are required in the parity odd 
sector of the model, while SM particles are in the parity even sector. Notably if 
this parity is exact, the lightest parity odd particle is stable, therefore the picture 
is consistent with the existence of the dark matter in our Universe. Finally, the 
production of the top partners and their decays into the top quark and an invisible 
particle are predicted as the signal of "Beyond the Standard Model" at the LHC. 

For the case of the MSSM, the pair production of the scalar top is not detectable, 
because the production cross section is too small. Supersymmetry is instead detected 
by the production of gluinos and squarks in the first generation, and the scalar top 
may appear as the decay product of the gluino [T2]. On the other hand, the cross 
section of the top partner pair production in the littlest Higgs model with T-parity 
is about 10 times larger than that of the scalar top, so that they may be detectable 
at the LHC [13J. The purpose of this paper is to provide a realistic estimate for the 
detection of the fermionic top partner at the LHC. 

This paper is organized as follows. In SecJU we review the littlest Higgs model 
with T-parity focusing on the top sector. In Secj3j we summarize the electroweak 
and dark matter constraints on the model. We find that the lower limit of the top 
partner mass is about 600 GeV, and the lightest T-odd particle is always much 
lighter than the top partner. The signal at the LHC is therefore two top quarks with 
significant transverse momentum and missing energy. 

In SecJHand[5l we discuss the top partner signature at the LHC. Reconstructing 
a top quark in the event is essential to identify the top partner. We describe the 
method to reconstruct the top quark in SecJH We apply a hemisphere analysis to 
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the signal reconstruction, and study the reconstruction efficiency of the top quark 
for both signal and ti background. In SecJSl we discuss the basic cuts to reduce the 
top quark background. The highest S/N ratio is obtained in the region where the 
effective mass is around twice the top partner mass. We also show the numerical 
results of our simulation study and find that the top partner signature would be 
significant over the background if the mass of the top partner is less than 1 TeV. 
Secj6] is devoted to discussion and comments. 

2 The Littlest Higgs Model with T-parity 

In this section, we briefly review the littlest Higgs model with T-parity focusing on 
the top sector of the model. The constraints on the model from WMAP observations 
and electroweak precision measurements will be discussed in the next section. For 
general reviews of the little Higgs models and their phenomenological aspects, see 

Refe.ua us]. 

2.1 Gauge-Higgs sector 

The littlest Higgs model [5] is based on a non-linear sigma model describing an 
SU(5)/SO(5) symmetry breaking. The non-linear sigma field, S, is given as 

S = e 2tn ^S , (1) 

where / is the vacuum expectation value associated with the breaking and is expected 
to be 0(1) TeV. The Nambu-Goldstone (NG) boson matrix n and the direction of 
the breaking S are 
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where we omit would-be NG fields in the n matrix. An [SU(2)xU(l)] 2 subgroup in 
the SU(5) global symmetry is gauged, which is broken down to the diagonal sub- 
group identified with the SM gauge group, SU(2)lxU(1)y. Due to the presence of 
the gauge interactions (and Yukawa interactions introduced in the next subsection), 
the SU(5) global symmetry is not exact, and particles in the n matrix become pseudo 
NG bosons. Fourteen (= 24 — 10) NG bosons are decomposed into representations 
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under the electroweak gauge group as lo © 3o © 2±i/ 2 © 3±i. The first two represen- 
tations are real, and become longitudinal components of heavy gauge bosons when 
the [SU(2)xU(l)] 2 is broken down to the SM gauge group. The other scalar s in the 
representations 2 ±1 / 2 and 3±i are the complex doublet identified with the SM Higgs 
field (if in Eq.(f2])) and a complex triplet Higgs field ($ in Eq.(f5])), respectively. The 
kinetic term of the E field is given as 



f 2 

£s = ^Tr 



f^E (f^E) f 



where 



f^E = c^E - i £ [^(WjE + EWj) + ^(B,E + SBj) 

3=1 



(3) 



(4) 



Here, W 3 - = WfQ* (Bj = BjYj) is the SU(2) j (U(l) j ) gauge field and gj{g\) is the 
corresponding gauge coupling constant. The generators Qj and Yj, are 



Qi 



( o a \ 


v j 

( \ 

, Y 2 = diag(2,2,2,-3,-3)/10 , (5) 
v o a * j 

where a a is the Pauli matrix. 

In terms of the above fields, the symmetry under T-parity [8] is defined as the 
invariance of the Lagrangian under the transformation: 



Yl = diag(3,3,-2,-2,-2)/10 



n «-> -QUQ (or E E = EnfiE^Er 



(6) 



where f2 = diag(l, 1,— 1,1,1). As a result of the symmetry, the gauge coupling gi 
(g[) must be equal to g 2 (g' 2 ), namely g x = g 2 = V2g (g[ = g 2 = V^g'), where g (g') 
is nothing but the coupling constant of the SM SU(2) L (U(l)y) gauge symmetry. 
The Higgs potential is generated radiatively P, [9| 



V(H, $) = A/ 2 Tr 



$t$ 



/i 2 fftif+^(iftif) 2 + ... . 



(7) 



The main contributions to /i 2 come from the logarithmic divergent corrections at 
1-loop level and quadratic divergent corrections at 2-loop level. As a result, /i 2 is 
expected to be smaller than f 2 . The triplet Higgs mass term, on the other hand, 
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receives quadratic divergent corrections at 1-loop level, and therefore is proportional 
to f 2 . The quartic coupling A is determined by the 1-loop effective potential from 
the gauge and top sectors. Since both /x and A depend on the parameters at the 
cutoff scale, we treat them as free parameters in this paper. 

Next, we discuss the mass spectrum of the gauge and Higgs bosons. This model 
contains four kinds of gauge fields, W®, W 2 , B 1 and B 2 , in the electroweak gauge 
sector. The linear combinations W a = (Wf + W£)/V2 and B = (B x + B 2 )/V2 
correspond to the SM gauge bosons for the SU(2)^ and U(l)y symmetries. The other 
linear combinations, W H = (W? - W£)/V2 and B H = (B 1 - B 2 ) /y/2, are additional 
gauge bosons, which acquire masses of 0(f) through the SU(5)/SO(5) symmetry 
breaking. After the electroweak symmetry breaking, the neutral components of W H 
and Bh are mixed with each other, and form mass eigenstates Ah and Zh- The 
masses of the heavy bosons are 

rriA H ^ OMg'f , m Z[I ^ m WlI ~ gf . (8) 

The mixing angle Oh between W H and B H is considerably suppressed. It is given 
by tan Oh — —g' v2 /{^gf 2 ), where v (~ 246 GeV) is the vacuum expectation value 
of the Higgs field. Thus the dominant component of Ah is Bh- Finally, the mass of 
the triplet Higgs boson $ is given by m% = Xf 2 = 2m 2 l f 2 /v 2 , where rrih is the SM 
Higgs boson mass. 

Under T-parity, the new heavy gauge bosons and the triplet Higgs boson behave 
as T-odd particles, while SM particles are T-even. As shown in Eq.flHD, the heavy 
photon is considerably lighter than the other T-odd particles. Stability of Ah is 
guaranteed by T-parity conservation, and it becomes a candidate for dark matter. 

2.2 Top sector 

To implement T-parity, two SU(2) doublets, q\ and q 2 , are introduced for each SM 
fermion doublet. Furthermore, two vector-like singlet top partners, U\ and U 2 , are 
also introduced in the top sector in order to cancel large radiative corrections to the 
Higgs mass term. Since we are interested in top partner production at the LHC, 
only the top sector is discussed here. For the other matter sectors, see Refs. [H [TB] . 

The quantum numbers of the particles in the top sector under the [SU(2) x U(l)] 2 
gauge symmetry are shown in Table [TJ All fields in the table are triplets under the 
SM SU(3) C (color) symmetry. Using these fields, the Yukawa interaction terms which 



5 



Qi 


(2, 1/30; 1,2/15) 


<?2 


(1,2/15; 2, 1/30) 


u L1 


(1,8/15; 1,2/15) 


u L2 


(1,2/15; 1,8/15) 




(1,8/15; 1,2/15) 


U R2 


(1,2/15; 1,8/15) 
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(1,1/3; 1,1/3) 







Table 1: The [SU(2)x U(l)] 2 charges for particles in the top sector. 



are invariant under T-parity and gauge symmetries turn out to be 



'-^-ijk^-xy 



2\/2 

-A 2 / (UliU r1 + U L2 U R2 ) + h.c. 



ky 



UR 



(9) 



where 



Qi 



( 9* \ 
V o J 



qi = -o 



ULi 



(10) 



The indices i,j,k run from 1 to 3 whereas x, y — 4,5. The coupling constant Ai is 
introduced as the top Yukawa coupling, while A2/ gives the vector-like mass term for 
the singlet fields. Under T-parity, qi and C7, transform as q\ <-> — q 2 and U\ —U 2 - 
Therefore, the T-parity eigenstates are given by 



1 1 

q± = Tq2) , U L± = -j=(U li tU L2 ) 



Ur± = ^=(U r1 tU R2 ) . (11) 



In terms of these eigenstates, the mass terms for these quarks are written as 
follows, 



-Al fU L+ + VU L+ Ur- X 2 f [U L+ U R+ + U L -Ur-) + h.c 



(12) 



The remaining T-odd fermion, q_, acquires mass by introducing an additional SO (5) 
multiplet transforming nonlinearly under the SU(5) symmetry. Therefore, the mass 
term of the g_ quark does not depend on Ai and X 2 . In this paper, we assume 
that the g_ quark is heavy compared to other top partners, and do not consider its 
production at the LHC For the g_ quark phenomenology, see Ref.[16j. 
The T-even states u + and U + form the following mass eigenstates 

ti, = cos (3 ul+ — sin (3 Ul+ , Tl = sin j3 ul+ + cos (3 Ul+ , 

tR = cos a ur — sin a Ur + , Tr = sin a ur + cos a Ur + , (13) 
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where sin a ~ Ai/(Af + A 2 ,) 1 / 2 , and sin/? ~ A 2 w/[(A 2 + A 2 .)/]. The t quark is identified 
with the SM top quark, and T is its T-even heavy partner. On the other hand, the 
T-odd fermions Ul- and Ur- form a Dirac fermion, T_. The masses of these quarks 
are given by 



It is worth noting that the T-odd states do not participate in the cancellation of 
quadratic divergent corrections to the Higgs mass term. The cancellation is achieved 
only by loop diagrams involving t and T quarks. 



In the previous section, four parameters are introduced in the gauge-Higgs and top 
sectors in addition to the gauge coupling constants (g and g') and the vacuum expec- 
tation value of the Higgs field (v). Those are m/,, /, X%, and A 2 . Since the top quark 
mass is determined by the combination of v, Ai and A2, the number of undetermined 
parameters is three. These parameters can be expressed by nth, rriAn-> and m^. 

In this section, we see that and rriA H are directly related each other thanks 
to the precise data of the WMAP observations. Therefore, it is possible to write 
down the Higgs mass as a function of the dark matter mass, = m^m^). The 
electroweak precision measurements provide a further constraint on rrih and mr_- 
They gives a lower bound on and a large mass difference between T_ and Ah- 

3.1 Constraint from WMAP observation 

First, we consider the WMAP constraint on the littlest Higgs model with T-parity. 
The dark matter, Ah, annihilates mainly into weak gauge bosons, W + W~ and ZZ 
through diagrams in which the Higgs boson propagates in the s-channel. Once we 
calculate the annihilation cross section, the thermal relic abundance of the dark 
matter is obtained by solving the Boltzmann equation. For the detailed calculation 
of the abundance in this model, see Refs.[9]-[TT]. To good accuracy, the dark matter 
abundance can be written as 






(14) 



3 WMAP and EW precision constraints 
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Figure 1: Contour plot of the thermal relic abundance of the dark matter, f^DM^ 2 , in the 
(jn,A H ,Tnh) plane. The thin shaded area is the allowed region from the WMAP observation 
at 2cr level, 0.094 < n DM h 2 < 0.129. 

where c is the speed of light, and v TB \ is the relative velocity between initial dark 
matters. Since the product of the cross section and the relative velocity, <7t> re i, and 
hence the dark matter abundance, Odm'* 2 , depend only on mA H and nth, the WMAP 
observation, Q-ouh 2 — 0.111, gives a relation between these parameters. 

In Figfr], the thermal relic abundance of dark matter is depicted as a contour 
plot in the (mA H ,rrih) plane. The thin shaded area is the allowed region from the 
WMAP observation at 2a level, 0.094 < tt DM h 2 < 0.129 [Hj. In the figure, there 
are two branches: the upper branch (U-branch) and the lower branch (L-branch). 
These branches are sometimes called "Low" and "High" regions, respectively |llj . 
In the U-branch, the Higgs boson mass is larger than twice the dark matter mass, 
rrih > 2m a h , while m^ < 2m a h in the L-branch. The Higgs mass is precisely 
determined by the dark matter mass up to a two-fold ambiguity by imposing the 
WMAP constraint. 

3.2 Constraint from electroweak precision measurement 

Next, we discuss the constraint from electroweak precision measurements. New 
physics contributions to electroweak observables come from radiative corrections, 
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U-branch L-branch 




400 600 800 1000 1200 1400 400 600 800 1000 1200 1400 

m T (GeV) m T (GeV) 

Figure 2: Constraints on mx_ and rriA H at 68% and 99% confidence level. At each point 
in these figures, the Higgs mass is determined to satisfy the WMAP constraint on the 
U-branch (left figure) and L-branch (right figure). 



since there is no tree-level effect due to T-parity. The constraint is sensitive to the 
masses of the Higgs boson and the top partner. Since the Higgs mass is directly 
related to the dark matter mass through the WMAP constraint, the electro weak 
precision constraint on rrih and itlt^ is translated into one on rriA H and . 

In order to obtain the constraint, we follow the procedure in Ref.jTS] using the S, 
T and U parameters [Tj5]. In that paper, it has been shown that main contributions 
to the parameters come from the top-sector and the custodial-symmetry violating 
effect from heavy gauge boson loops. However, in Ref.[10j. it has been shown that 
the latter contribution is negligibly small compared to the former one. Therefore, 
we consider only the top sector contribution in order to obtain the constraint. For 
the detailed expression of the top sector contributions, see Ref.|18j. 

In FigfJl constraints on and rriA H at 68% and 99% confidence level are 
depicted. At each point in these figures, the Higgs mass is determined to satisfy 
the WMAP constraint on the U-branch (left figure) and L-branch (right figure). To 
obtain the constraint, we have used three experimental values: the W boson mass 
(m w = 80.412 ± 0.042 GeV) and the weak mixing angle (sin 2 ^ e ff pt = 0.23153 ± 
0.00016) [20], and the leptonic width of the Z boson (Tj = 83.985 ± 0.086 MeV) 
[2"T| . We have also used the fine structure constant at the Z pole (a _1 (m^) = 
128.950 ± 0.048) and the top quark mass (m t = 172.7 ± 2.9 GeV) [22]. 

As seen in these figures, the lower bound of the T_ quark mass is about 600 GeV, 
and the mass difference between rriT_ and rriA H is larger than 500 GeV0. In the 

: As shown in the previous subsection, the dominant annihilation mode of dark matter for the 
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excluded region painted black, there is no combination of the parameters Ai, A2 and 
/ which can give the correct top quark mass via Eq. (|14p . The parameter region with 
large rriA H corresponds to the region with heavy rrih due to the WMAP constraint. 
The heavy Higgs contribution to the T-parameter can be cancelled by those from 
the top partner, T, if its mass is tuned appropriately. When rriA H ~ m h /2 increases, 
the cancellation can be achieved only for small region of as we can see in these 
figures. 

4 T_ quark production at the LHC 

In this section, we study the signature of T_ quark production at the LHC. Following 
the discussion in the previous section, we first discuss properties of the T_ quark and 
present a few representative points used in our simulation study. We next consider 
a top reconstruction for the signal using a hemisphere analysis. Finally, we address 
the background to this process, which comes from SM ti production. 

4.1 Properties of the T_ quark 

Due to the T-parity conservation, the T_ quark would be pair produced from proton- 
proton (pp) collisions at the LHC. Since the T_ quark has a color charge, it is 
produced dominantly through SU(3) C interactions at the LHC. The production cross 
section depends only on its mass, my_. Unlike the scalar top in the MSSM, the 
T_ quark is a Dirac fermion. Hence its production cross section at the LHC ranges 
between 0.1-1 pb when m*r_ is less than 1 TeV as shown in Figj3] (left figure). 

The decay process of the T_ quark is simple, because only Ah and Zh (Wh) are 
T-odd particles lighter than the T_ quark. Furthermore, the T_ quark is almost an 
SU(2)x singlet. The interactions relevant to the decay are 

£ = i-^- cos 9hT-A h (sin /3 Pt, + sin ckPr) t 
2g' 

+ i— sin 6 H T_ft H (sin j3P L + sin o;Pr) t . (16) 

As seen in equation ([TBI) , the decay mode T_ — > Z H t is highly suppressed by sin 2 6 H , 
therefore T_ decays dominantly into the stable Ah and the top quark. In Figj3] (right 
figure), the decay width is depicted as a contour plot in the (m T _,mA H ) plane. It 
shows that the width is typically several GeV. 
relic abundance is Ah Ah — > W + W~ , implying that rriA H > mw- 
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Figure 3: Production cross section of the T_ quark at the LHC as a function of (left 
figure). Contour plot of the T_ decay width in the (mT_,mA H ) plane (right figure). 





m T _ (GeV) 


m AlI (GeV) 


a(pp^T_T.+X) (pb) 


I 


600 


100 


0.940 


II 


700 


125 


0.382 


III 


800 


150 


0.171 


IV 


900 


175 


0.0822 



Table 2: The model points for our simulation study. The production cross section for the 
T_ quark is also shown. 



The signal of T_ quark production is a top pair (tt) with significant missing 
transverse momentum. In this paper, we assume that the branching ratio of the 
process T_ — > A^t is 100%. For simplicity we also assume that the other extra 
matter fermions do not contribute to the signal. The model points we choose for the 
simulation study are listed in Table [2j The production cross sections are also shown 
in the table, which were obtained by the CompHep code |23j using the "CTEQ6L1" 
parton distribution function [21] and the scale of the QCD coupling set to be Q = 
m T _§. 

In order to generate parton level events, we calculated the process pp — > UAhAh 
directly using the CompHep code, keeping diagrams relevant to the process through 

2 The cross sections are a factor of 2 smaller than those listed in 13J. We found their cross 
section is obtained at the scale of mz- 
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on-shell T_T_ production. We generated 100,000 events for each model point to 
study their distributions. The parton level events were interfaced to HERWIG [22] 
for fragmentation, initial and final state radiations, and hadronization. The effect 
of the top polarization is not included in our simulation. The detector effects were 
simulated by the AcerDET code [22] • This code provides a simple detector simulation 
and jet reconstruction using a simple cone algorithm. It also identifies isolated 
leptons and photons, finds b and r jets, and calculates the missing momentum of the 
events using calorimeter information. 

Before going to the discussion of top reconstruction, we define two important 
quantities which are frequently used for new physics searches at the LHC One is 
the missing transverse energy, -Ermiss = {Ptx + -Pt J /) 1 ^ 2 > "which is important for signal 
of the models with a stable dark matter candidate such as the MSSM or the little 
Higgs model with T-parity. Here Pt% is the sum of transverse momenta measured 
by the calorimeter. The other is the effective transverse mass, 

M eS =J2PT+ E PT+ E PT + ^Tmiss, (17) 

jets isolated leptons isolated photons 

where we require that the pseudo rapidity is less than three (rj < 3) and the trans- 
verse momentum is pt > 50 (10) GeV for each jet (lepton/photon). The effective 
transverse mass is a good quantity to measure the mass scale of a produced particle. 

4.2 Top reconstruction 

Now we move on to the discussion of top quark reconstruction for the signal. The top 
quarks produced from T_ quarks have significant p?- Since the T_ quarks must be 
produced in a pair, we expect two separate jet systems originating from the two top 
quarks. We therefore use the hemisphere analysis [27] for the event reconstruction. 
Two hemispheres are defined in each event, and high p? jets, leptons and photons 
are assigned to one of the hemispheres. Specifically, 

• Each hemisphere is defined by an axis Pi (i = 1,2), which is the sum of the 
momenta of high p? objects (jets/leptons/photons) belonging to the hemisphere 
i. Only the jets with pt > 50 GeV and leptons/photons with pt > 10 GeV 
are assigned to the hemisphere in order to reduce the contamination of QCD 
activity. 

• High pt objects k belonging to the hemisphere % satisfy the following condition 

d(p k ,Pi) <d(p k ,P 3 ) , (18) 
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Figure 4: Distribution of a two jet invariant mass, m max (jj) (left figure), and that of 
m m i n (jjj) (right figure) in the pp — > T_T_ — > ttAjjAu process for point III 



where i and j are the indices of the hemispheres, and the function d is defined 

as 

d( Pk ,Pi) = (Ei - \Pi\ cos9 ik ) , (19) 

{Ei + E k y 

where 9i k is the angle between Pi and p k . 

To find the axis satisfying the above conditions, we take following steps. (1) 
We take the highest px object i (jets/leptons/photons) with momentum p iy and the 
object j with largest AR\pj\, where AR = [(A<p(i, j)) 2 + (Arj(i, j)) 2 ] 1 / 2 . We take 
Pi and pj to be the seeds of the hemisphere axes, namely, P[ n = pi, P™ = pj. (2) 
Each object with momentum p k is assigned to the hemisphere i, if d(p k , Pl n ) < 
d(p k ,Pj n ). (3) We then define new P- n (i = 1,2) as the sum of the momenta of 
the objects belonging to the hemisphere i. (4) We repeat the processes (2) and (3) 
until assignment converges. In this paper, we denote the hemisphere seeded from 
the highest px object as "hemisphere 1" and the other as "hemisphere 2". 

If there are more than two jets in a hemisphere, we can calculate the maximum of 
the two jet invariant masses for all combinations of jets in the hemisphere, m maiX (jj). 
In FigJH (left figure), we plot this for the hemisphere 1 using the point III in Table 
[2J The distribution has two peaks around m(jj) ~ 80 GeV and 130 GeV. The lower 
peak corresponds to the jet pair arising from the W decay, while the second peak 
comes from the combination of one of the two jets from a b quark and one of the 
partons from the W decay. 
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pp — > T-T- 


pp — > tt 


m min (jjj)i < 200 GeV or m min (jjj) 2 < 200 GeV 


22.9% 


15.3% 


m m in(jjj)i < 200 GeV and m min (jjj) 2 < 200 GeV 


1.6% 


0.17% 



Table 3: Fraction of T_T_ and ti events with m m i n (jjj) < 200 GeV in one of the hemi- 
spheres, and both of the hemispheres. The probability to reconstruct the top quark in 
both of the hemisphere is significantly small for the tt production. 



When a hemisphere contains more than three jets, we can also define the mini- 
mum three jet invariant mass in the hemisphere, rn min (jjj), where two of the three 
jets are those which give m max (jj). In FigJH (right figure), the distribution peaks 
at the input top quark mass 175 GeV, clearly showing that the hemisphere analysis 
can group the jets from the top quark correctly. In the following, we often require 
the "top mass" cut, m m i n (jjj) < 200 GeV for at least one of the hemispheres. 

The efficiency to find at least one top quark candidate in a signal event is mod- 
erate, about 20%. In FigJU we have a long tail at m min (jjj) > 200 GeV, which 
consists of events with additional jets or miss- assignment of jets. If we optimize 
the top search strategy after the hemisphere reconstruction, we may increase the top 
reconstruction efficiency, but we do not study this possibility in this paper. It should 
be noted that the reconstructed m m i n (jjj) distribution in the hemisphere analysis is 
not biased, because we do not assume the existence of the top quark in the event. 



4.3 Background 

The most important background comes from the pp — > ti+X process in the SM. The 
tree level production cross section of the top quark is 400 pqj. We have generated 
events corresponding to 50 fb _1 for this study. The cross section after very weak 
cuts, M eff > 400 GeV, E Tmiss > 100 GeV, n ]ets (p T > 100 GeV) = n 100 > 1 and 
njctsipr > 50 GeV) = n 50 > 2 is about 16 pb, which is still higher than the signal 
cross section by more than a factor of 10. 

The efficiency to find the top quark candidate in a hemisphere becomes lower 
for the tt production process if large -Exmiss is required. This is because the W 
boson from the top quark decay must decay leptonically to give such a high -Exmiss 
to the event. This can be seen in Table [31 where we have listed the probability 
to find three jets with m min (jjj) < 200 GeV in one, or both, of the hemispheres. 



3 The NLO cross section is 800 pb 
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In FigJHJ we present the m min (jjj) distribution in a two dimensional plot, where 
the x (y) axis corresponds to m m i n (jjj) in hemisphere 1 (hemisphere 2). While the 
distribution peaks around m t in both of the hemispheres for the signal (left figure), 
the distribution for ti production scatters over the plot (right figure). This means 
that at least one of the top quarks has to decay leptonically, and therefore some of 
the jets must be additional QCD jets in ti production. 

If we require a top candidate in both hemispheres, the ti background reduces to 
the 30 fb level, however the reconstruction efficiency for the signal also decreases 
due to additional jets in the final state and the overlap of jets. In the next section, 
we look for the excess of signal events over the ti background for the events with 
m m i n (jjj) < 200 GeV in at least one of the hemispheres, for the various different 
sample points. In this case, the ti background cross section is around 16 pb x 0.153 
~ 2.5 pb, which is many orders of magnitudes larger than any other irreducible 
backgrounds listed in [T5] , 

In Ref . |29] . it has been shown that the quark production is not the dominant 
background in the inclusive study of SUSY processes for the 0-lepton + -Exmiss chan- 
nel. The backgrounds from W and Z boson production are about as important as 
that due to ti production. However, one should be able to reduce these backgrounds 
significantly by requiring a reconstructed top and b tagged jets. We therefore do not 
study these in this paper. 

5 Discovery of the T_ quark at the LHC 

In this section, we investigate the possibility to find the T_ quark at the LHC. First, 
we discuss the separation of the signal from the background using their different 
kinematic properties. We find the region in the -&rmiss and M e g plane where the 
signal best dominates the background. Next, we calculate the statistical significance 
for T_ quark discovery at the LHC, and show that the significance exceeds seven 
for all sample points. We also calculate the M^ 2 variable [30J to investigate the 
possibility to extract information about mr_ from the signal. Finally, we comment 
on the differences in production and decay distributions between T_ and scalar top 
signals. 
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Fi gure 5: Reconstructed ?72 mm (jjj) distributions for pp — ► T—T— for point III (left figure) 
and pp — > ti (right figure) processes. We take i^Tmiss > 400 GeV and nioo > 1 and 
^50 > 2 for the ti events. The efficiency to find the top candidate in each hemisphere 
(wmin (jjj)i < 200 GeV) is low in the ti production compared to the signal process. 



5.1 Separation of signals from backgrounds 

In order to separate T_T_ signals from ti backgrounds, we further study E^ miss 
distributions for a given (high) M e fj interval. In Fig® we show the signal and 
background distributions in the M e g and Ermiss plane. The Ermiss distribution of the 
signal peaks near its maximum (~ 0.5M e g) for significantly large M e g. This feature 
is common in processes where new particles are pair produced and each decays into 
a stable neutral particle and other visible particles [2TJ. This can be understood as 
follows. Since new heavy particles are produced mostly near the threshold at the 
LHC, the velocity of the T_ quark is low in the transverse direction. Then, -Ermiss 
and M e ff are maximized when the decay configuration of the two T_ quarks is such 
that the two top quarks from the decay go in the same direction in the rest frame 
of the T_T_ system, as illustrated in FigJTl In this case, M eS ~ 2£^ Tmiss ~ 2my_. 

On the other hand, the kinematics is totally different for the background ti 
distribution. As shown in FigJTl neutrinos arising from top decays are collinear to 
the direction of the parent top quark, if the center of mass energy in the collision is 
much higher than the top quark mass. Therefore, E^ miss ~ 0.5M e fj is kinematically 
disfavored. This can be seen in Figj6] (right figure). As M e g increases, the fraction 
ErmissA^eff is significantly reduced. 

Using the nature of these production and decay kinematics, we can find the 
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Meff(GeV) Meff(GeV) 

Fi gure 6: -E/rmiss versus M e g distributions for the T_T_ production at the LHC for point 
III (left figure), and the tt background (right figure). Normalizations in both figures are 
arbitrary. 

kinematical region with good separation between signal and background. We restrict 
M eff to a certain large value, so that we see the bump of the signal event in the Ermiss 
distribution. In FigJHl we show the signal and background E Tmiss distributions for 
points I to IV (top four figures). Each plot corresponds to the integrated luminosity 
JdtC = 50 fb _1 . We restrict M cS to 2m T _- 200 GeV < M eS < 2m T _ for m T _ = 600 
and 700 GeV, and 2m T _- 300 GeV < M eS < 2m T _ for m T _ = 800 and 900 GeV, so 
that -ETmiss becomes maximal and the signal rate is still reasonably high. Here the 
top mass cut, m min (jjj) < 200 GeV, is required for at least one of the hemispheres, 
which reduces the background by factor of 5 and the signal by factor of 3 compared 
to the case where no cut is applied on m m i n (jjj). In the bottom four figures in 
FigJHl we further require that the events have no isolated leptons. The isolated 
lepton is produced by the leptonic decay of the top quark. The it background is 
reduced by a factor of two by this cut, with no significant reduction of signal events. 
Each distribution shows the clear excess of events over the (exponentially decreasing) 
background. 

The excess is less prominent for smaller M e g bins. We show the Ermiss distribution 
for point II with 1000 GeV < M eS <1200 GeV in FigJHl (left two figures). While 
the number of the signal increases by factor of two, the it background increases 
more than factor of three. On the other hand, if one increases M e g, the number of 
the signal reduces rather quickly, and we find that the events in this region do not 
contribute much to the discovery of the T_ quark. See FigJHl (right figure) for the 
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Figure 7: T_T_ production which gives the largest E^ B 
ti production with neutrino emissions. 
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Signal/BG 


Signal/BG 


Signal/BG 


(GeV) 


(GeV) 


(GeV) 


(0-lepton with top cut) 


(with top cut) 




600 


1000 


400 


842/106 


1053/313 


3336/1304 


700 


1200 


450 


263/54 


332/114 


1284/582 


800 


1300 


500 


208/28 


249/57 


874/417 


900 


1500 


550 


93/7 


105/16 


397/203 



Table 4: The signal to background ratio at the sample points I to IV. We take a region 
with Mtff n < M e fj < 2rriT_ and £/rmiss > -E'x'miss- The ra ti° is best if one requires the top 
cut and vetoes isolated leptons. 



distribution with 1400 GeV < M c g <1600 GeV. Note that the events in the region 
M eS 3> 2m T _ arise from highly boosted T_ quarks, therefore the fraction £ , Tmiss /M eff 
decreases, making the signal distribution less prominent over the background. 

The numbers of the signal and background events in the signal region are shown 
in Table HI Here we take the same signal region as that of FigJHl The signal to 
background ratio with the top cut is more than 3 for all sample points. Thus, it is 
clearly shown that the signal dominates in the region where -Epmiss ~ 0.5M e fj. 



5.2 Statistical significance for the T_ quark discovery 

As shown in FigJHl the ti background reduces rather quickly at large i?Tmiss- Near 
the end of the distribution of -Ermiss; the distribution is dominated by the signal in 
Fig. [SJ On the other hand, the signal distribution decreases quickly after its peak. 
We therefore fit the decrease of total distribution near the maximum value ~ 0.5M e fj 
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Figure 8: i^Tmiss distributions for points I to IV from left to right, with 2m^_ — 200 
(300) GeV < M cS < 2m T _ for m T _ < 700 (> 800) GeV. The red and blue histograms 
are for the signal and background distributions, and the green histograms are the sum of 
these. m m i n (jjj) < 200 GeV is required at least in one of the hemispheres. In the bottom 
four figures, only the events without isolated leptons have been used. 
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Figure 9: £rmiss distributions for point II 
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with 1000 GeV < M eff < 1200 GeV (left two 
(right two figures). rn m i n (jjj) < 200 GeV is 
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m T _ (GeV) 


h 


Ah 


a (GeV) 


Nun 


AE hin (GeV) 

Dill \ vj t y 


h/Ah 


600 


190 


14.7 


38.5 


7 


20.8 


13 


700 


160 


6.8 


43 


6 


25.6 


24 


800 


46 


4.9 


66 


6 


30.4 


9.4 


900 


17 


2.4 


76 


6 


35.2 


7.1 



Table 5: Fit of the -&rmiss distribution to the Gaussian function near 0.5M e fj. Here, N^- m 
is the number of bins used for the fit and AE^in is the bin size. 



to the Gaussian function 

-^(-^Tmiss) 



hexp 



(3 



Tmiss 



E, 



( av g) \2* 
Tmiss J 



2(X 2 



(20) 



The result is summarized in Table [5j where we have used six or seven bins with the 
bin size (A£ bin ) between 20.8 and 35.2 GeV. These fits give A X 2 /{N hin - 3) ~ 1. 

The statistical significance of the signal is given by h/Ah, where h is the height 
of the Gaussian distribution and Ah is its error. The significance is more than seven 
for all sample points. When we fit the background distribution to the same Gaussian 
function, we obtain a = a t t as 102(125) GeV for the fit above -Ei-miss > 304(328) GeV 
for 1000 GeV < M cS < 1200 GeV (1200 GeV < M cS < 1400 GeV) respectively. The 
X 2 of the fit is Ax 2 /(iV bin - 3) = 1.1 (1.2). We find a < a th therefore, the detected 
edge is clearly inconsistent with the ti background distribution. 

We found that the signal is prominent over the background when M e g- is restricted 
to the region near 2m^_ . We therefore calculate the Mt2 variable for the events near 
the bump. We use "Cambridge M T2 " , which is defined as 



T2 



mm 



maxK_ (pf , p? H ),m 2 T _ (pf, p^)} 



(21) 



It is a function of the transverse momenta and masses of the two visible particles, 
-Prmiss, and the mass of the invisible particle. Mt2 is sensitive to the mass difference 
rriT_ — rriA H , but is not sensitive to the overall mass scale. We take events where 
there are more than three jets with >50 GeV in each hemisphere. We calculate 
Mt2 by taking the visible particle momenta as the sum of the three jet momenta 
used to calculate the minimum three jet mass m m i n (jjj) of the hemispheres, and we 
fix fiiAj = 150 GeV. We also require m mm (jjj) < 200 GeV for at least one of the 
hemispheres. The distributions are shown in FigflQl We found a positive correlation 
to the mass of the top partner, however the number of events that can be used for 
this analysis is small. 
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m T _= 600 GeV m T _= 700 GeV m T = 800 GeV 
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Figure 10: Distribution of the Mtz variable for the events in the signal region defined in 
Table H 

m rain(jjj) < 200 GeV is required in at least one of the hemispheres, nvr_ — 600, 
700, 800 GeV from left to right. 

5.3 Difference between T and scalar top signals 

We discuss the differences in the production and decay distributions between the 
T_ quark and scalar top (t) signals. It is impossible to find the process pp — > tt* 
followed by the decay t — ► txl at the LHC, because the production cross section is a 
factor 10 smaller than the T_ quark production cross section. However, we can still 
learn something from the comparison. 

Here we take = m 4 ~ = 800 GeV and rriA H = m^o = 150 GeV. We took low 
energy parameters of the MSSM so that all particles except stop and the LSP (x?) 
are too heavy to be accessible at the LHC, and t decays into t and Xi with 100% 
branching ratio. The events for the pp — > tt* — ► ttXiXi process are generated by the 
CompHep code and interfaced to HERWIG. 

In Figure [TT] (left figure), we show the pt distributions of T_ and t production 
at the LHC. The cross section is dominated by the low pt component for the T_ 
quark. On the other hand, the t production cross section is dominated by P wave, 
therefore the px distribution is broadened. The production cross sections peak at 
350 and 600 GeV for the T_ quark and t, respectively, where (3 ~ 0.4 and 0.6. Both 
cross sections are kinematically suppressed near the peak by (3. Especially in the 
stop case, the squared amplitude is proportional to f3 2 due to the P wave production. 
The cross section beyond the peak is suppressed by the quickly decreasing quark and 
gluon parton distribution functions. If this pt distribution is reconstructed from the 
tt distribution, this indirectly suggests that the produced particle is a fermion. 
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Figure 11: The px distribution of T_ (solid line) and t (dashed line) at the LHC (left 
figure). The 10 5 events are generated by the CompHep code. We set rriT_ = mi = 
800 GeV. The E/rmiss distribution for T_ (solid line) and t (dashed line) productions with 
1300 GeV < M c fj < 1600 GeV (center and right figures). For the right figure, we required 
m m i n (jjj) < 200 GeV for at least one of the hemispheres. 



The difference of the distributions affects the -Ermiss distribution for fixed M e ff as 
shown in Fig {11] (center and right figures). Here, the center figure shows the event 
distribution with 1300 GeV < M cS < 1600 GeV, and, in the right figure, when the 
top cut is required for at least one of the hemispheres. The -Ermiss distribution of the 
T_ events has a sharper peak near 0.5M e ff, while the distribution is broader for t. 
This is because t is produced with higher p?, therefore the top quark has, on average, 
a larger energy. Therefore the average fraction E^ miss /M e ff becomes smaller. 

Finally we comment on the top reconstruction for the supersymmetric signature. 
In supersymmetric models, the top quarks mostly arise from gluino decays through 
processes, g — > tt, g — > bb — > tbx + ■ The top quark in the event is hardly seen in the 
m min (jjj) distribution due to the additional high p T jets and leptons. In previous 
studies, the top quark is therefore searched for by looking for the jets consistent 
with the top decay kinematics, namely m(jj) ~ 80 GeV and m(bjj) ~ 175 GeV by 
looking for the combination of jets i, j, and k, which minimizes A% 2 defined as 

2 _ (m(i,j)-m w ) 2 (m(i,j,k) - m t f 
AX ~ (Am w y + (Am,) 2 ' 1 ' 

where Am^ and Am t are the expected errors on W and top mass reconstruction. 

The efficiency to find the top quark in the events becomes higher if we look 
for the jet combination consistent with top decay kinematics, however, in this case, 
we occasionally find top candidates in events which do not contain parton level 
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top quarks, increasing backgrounds from the other SM processes (32]. The signal 
distributions may also be distorted by such a procedure. This is particularly the case 
when the number of jets in the events is large. A scheme to subtract the accidental 
background from the top decay kinematical distributions has been developed |12j . 

6 Discussion and outlook 

In this paper, we have studied the phenomenology of the top partner in the littlest 
Higgs model with T-parity. This model predicts a stable neutral massive gauge boson 
{Ah), and a relatively light top partner (T_). The EW precision measurements and 
dark matter observations constrain parameters of the model, /, m^, Ai and A2. 
We find rriA H and are strongly constrained by the relic density constraint, if 
coannihilation processes with other T-odd particles are not efficient. Combined with 
the EW precision measurements, we have found the lower limit of m,T_ as a function 
of /, and clarified the fact that the top partner mass must be above 600 GeV and 
the mass difference between T_ and Ah is large. The signature of T_ quark pair 
production at the LHC is therefore high p? top quarks from T_ — > tAn and missing 
Pt coming from the Ah- We have studied the signal and background distributions 
at the LHC. 

Unbiased reconstruction of the top quark is essential to establish the existence of 
the top partner without increasing accidental backgrounds. We apply the hemisphere 
analysis, which is an algorithm to find two axes in events originating from particles 
produced in pairs. We have found that the appropriately defined three jet invariant 
mass in a hemisphere shows a clear peak at the top quark mass with a tail due to 
mis-reconstructed events. The three jet invariant mass is calculated from the pair 
of jets giving the maximum invariant mass in a hemisphere and a jet giving the 
minimum three jet invariant mass when combined with this jet pair. The efficiency 
to reconstruct at least one top quark, i.e. a hemisphere with m min (jjj) < 200 GeV 
in T_T_ events, is reasonable ~ 20%. 

The dominant background to the T_T_ process is ti production, whose production 
cross section with -Exmiss > 400 GeV is O(10) or more times larger than the signal 
cross section. We study the -Exmiss distribution as a function of M e g. Due to a 
simple kinematic reason, the separation of the signal from background becomes best 
when M e g ~ 2m^ . A pseudo-edge structure of the signal can be observed over the 
ti background in that region, and the significance of the signature turns out to be 
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larger than 7a for the case of m T _ < 900 GeV. 

The reconstructed signal is clearly different from that of SUSY events. In the 
case of SUSY, the production cross section of it* is rather small and is not detectable. 
Scalar top quark may be produced from gluino decay. In that case, additional jets 
are also produced, and the top cut m min (jjj) < 200 GeV for the jets in a hemisphere 
can hardly be satisfied. 

Results of our simulation given in this paper should be regarded as an order of 
magnitude estimate. We have used the simple smearing and jet finding algorithms 
provided by the AcerDET code. The E Tmiss and jet energy resolution could be dif- 
ferent at the ATLAS and the CMS detectors in the LHC environment. Note that 
the energy resolution assumed by the AcerDET is closer to the ATLAS detector ( 
The default is AE/E = 50%/ VE, while AE/E = 50(100)%/^ for the ATLAS 
(CMS) detectors, respectively). Furthermore, the background from tt production is 
estimated by the HERWIG code in this paper, where a(ti) = 400 pb. However, it is 
well known that the tt cross section receives large NLO corrections of around a factor 
of two. The multi-jet final state tt + jets should also give an important background 
to the events with -Ermiss- For the parameters we have studied, S/N is high in the 
signal region. Moreover, the signal has a peak structure in the E^ miss distribution. 
Therefore, we believe that our result is stable against additional sources of back- 
ground and a factor of two increase of the background. In addition, tt production 
with a multi-jet final state has a better chance to have a large three jet invariant 
mass in a hemisphere, because events have additional jets which are not kinemati- 
cally constrained by the top mass shell conditions. It is likely that the process tt + 
jets would be rejected by the top cut m m i n (jjj) < 200 GeV. It would be interesting 
to study the processes tt + jets and T_T_ + jets to check this expectation explicitly. 

In this paper, we do not study mass reconstruction in detail. The mass difference 
rriT_ — rriA H could be extracted from the M e s — Ermiss distribution if the background 
distribution can be calibrated precisely. We have also looked into the M T 2 distribu- 
tion for the sample with at least one reconstructed top quark, which shows positive 
correlation with the mass difference. The fit may be improved if the efficiency to 
reconstruct the top quark is improved by selecting the jets consistent with the top 
decay kinematics. The combinatorial backgrounds may be reduced if we select the 
top candidate among the jets in a hemisphere. 

The existence of the top partner and dark matter associated with the T-parity 
symmetry is an important feature of physics beyond the Standard Model as we dis- 
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cussed in the introduction. The collider phenomenology of the top partner therefore 
deserves more realistic studies. We hope more realistic studies will be performed by 
the ATLAS and the CMS groups in the future. 
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